Longitudinal Zeeman slowers based on permanent magnetic dipoles 
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Longitudinal Zeeman slowers composed of arrays of compact discrete neodymium magnets are 
proposed. The general properties of these slowers, as well as specific designs of short spin-flip 
Zeeman slowers for Sr and Rb atoms are described. The advantages of these slowers are their 
simplicity, low cost and absence of consumed electrical power and corresponding water cooling. The 
smoothness of the magnetic field together with ease of adjustability makes it possible to operate 
these slowers near the theoretical limits of deceleration, making them more compact and efficient. 

PACS numbers: 32.80.Lg, 32.80.Pj, 39.10.+J 



I. INTRODUCTION 

A Zeeman slower is a special magnet, which allows con- 
tinuous laser cooling of thermal atomic beams with just a 
single laser beam of fixed frequency. The purpose of this 
magnet is to compensate the Doppler frequency shift of 
the cooling atomic transition by the spatially-dependent 
Zeeman frequency shift, which provides a continuous res- 
onant interaction between the cooling laser light and the 
decelerated atoms for all atomic velocities within the cap- 
ture range of the slower. As the most intense sources 
of cold atoms, ovens combined with Zeeman slowers are 
widely used in different atomic physics experiments, in- 
cluding producing Bose Einstein Condensates p], [|[ and 
optical atomic clocks Q, which demand a large number 
of cold atoms. 

A standard longitudinal Zeeman slower [4| uses a spe- 
cial current-carrying coil to create a proper spatial dis- 
tribution of the magnetic field along its axis. 

In our previous work Q , Zeeman slowers based on ar- 
rays of permanent magnetic dipoles (MD) were proposed. 
The first of these spin-flip transverse Zeeman slower for 
Sr atoms @, based on compact neodymium permanent 
magnets, was recently implemented in the NPL Sr op- 
tical clock and has already demonstrated its efficiency, 
convenience and reliability. In spite of the simplicity of 
the transverse MD-Zeeman slowers (where the magnetic 
field is directed transversely to the direction of the atomic 
and cooling laser beams) such slowers have some limita- 
tions. First, the transverse Zeeman slower demands twice 
as much laser power, compared to standard longitudinal 
slowers. For linear polarization of the cooling light, only 
one of its two circular-polarization components of differ- 
ent helicity is in resonance with the corresponding Zee- 
man shifted transition of the slowed atoms. Second, the 
spin-flip (the slower, in which the magnetic field is chang- 
ing its sign) transverse Zeeman slower is working only for 
atoms without (or with small) Zeeman splittings of the 
ground internal state, like Sr or Yb, for example. For 
atoms with essential splitting of the ground state, like 
Li, Na, Rb, Cs, Cr the linear polarization of the cooling 
light of the transverse Zeeman slower prevents it from op- 
tically pumping atoms to the correct (slowed) magnetic 
substate in the region where the the magnetic field turns 



to zero, which leads to loss of the decelerated atoms at 
that region. This last fact has been confirmed by the 
recent tests of the transverse MD-Zeeman slower for Li 
and Cr atoms 0j- 

Recently a longitudinal Zeeman slower for Rb atoms, 
based on an 8-bars Halbach configuration of permanent 
magnets, has been successfully demonstrated |8[. 

This paper considers longitudinal spin-flip Zeeman 
slowers, which are using just two or four rows of standard 
cylindrical neodymium magnets, configured according to 
the original general proposal @ . The simplest of these Sr 
longitudinal slowers has already been successfully tested 
on the NPL optical clock setup @. 

In the first chapter of this paper, a short revision of 
the general theory of Zeeman slowers is given. The sec- 
ond chapter describes a general approach to the design of 
longitudinal MD-Zeeman slowers based on arrays of mag- 
netic dipoles (MD) as well as the two designs of slower 
for Sr atoms. The fourth chapter gives an example of a 
design of a spin-flip longitudinal slower for Rb atoms. Fi- 
nally, in the conclusion the main results are summarized. 



II. GENERAL THEORY OF A ZEEMAN 
SLOWER 

A Zeeman slower decelerates atoms by resonant scat- 
tering of photons from a counterpropagating laser beam. 
It is supposed that the monochromatic laser light of fre- 
quency ui resonantly interacts with some isolated and 
closed atomic transition of frequency ujq, such that the 
corresponding interaction can be described with a two- 
level model of an atom. In addition, it is assumed that 
the atom is placed in the magnetic field of a Zeeman 
magnet, which is characterized by a spatially dependent 
amplitude B(z). We assume also that this magnetic field 
leads to a Zeeman shift of the frequency of the atomic 
transition of fi'B(z)/H, where \i! is the magnetic moment 
for the atomic transition. The average spontaneous light 
pressure force on the atoms in the direction of propaga- 
tion of the laser beam, given by 
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FIG. 1: The solid line gives the dependence of the relative 
light force, characterized by the coefficient e (see the text), 
on the relative effective frequency detuning of the cooing light 
Aeff/T. The dashed line gives the same dependence of the 
derivative of the force. 



where k is the wave vector of light, sq(z) is the local on- 
resonance saturation parameter of the atomic transition, 
r is the linewidth of the transition, <5 = u — lo is the 
laser frequency detuning and v is the velocity of the atom. 
The velocity dependence of the force is determined by the 
effective frequency detuning A e ff — So + kv + n'B(z)/h, 
which includes the Doppler shift kv of the atomic fre- 
quency. The dependence of the relative magnitude of 
the force on the effective frequency detuning A e //, for 
sq(z) = 2, is shown in fig. 1. 

The maximum value of the local deceleration, provided 
by the force, is achieved at exact resonance, A e /^ = 0, 
and is given by 
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where m is the mass of the atom. Although the maximum 
deceleration can be used to estimate the shortest possi- 
ble length of a Zeeman slower, this can not be realized in 
practice. At exact resonance an equilibrium between the 
inertial force of the decelerated atoms and the light force 
is unstable and any slight increase of the atomic velocity 
due to imperfection of the magnetic field distribution or 
spontaneous heating of atoms will lead to a decrease of 
the decelerating force and subsequent loss of atoms from 
the deceleration process. In practice the deceleration of 
atoms is realized at a fraction of the maximum decelera- 
tion 
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where e < 1 ||. Note that the coefficient e in our case 
corresponds to the ratio between the reduced local ac- 
celeration and the maximum possible acceleration at the 
same location, which is a function of the local satura- 
tion parameter sq(z). In a standard treatment [ToL Ull] 
a similar coefficient -q appears, which relates the actual 
deceleration to the maximum possible deceleration at in- 
finite intensity of laser light. The two coefficients e and r\ 



are related as e = r\so{z)j(\ + s (zj). Therefore, the con- 
stant coefficient e corresponds to a variable coefficient r\, 
the value of which depends on the local intensity of light. 

When the local deceleration a(z) is less than a max , 
which corresponds to e < 1 the decelerated atoms stay 
on the low-velocity wing of the Lorentzian velocity pro- 
file of the light force (1) and undergo stable deceleration 
(see fig. I). The corresponding frequency detuning of the 
cooling light from the optical transition of the decelerated 
atoms at the active part of a Zeeman slower is determined 
by 



soO)) : 
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In a general case when the saturation parameter sq(z) 
changes due to convergence of the cooling laser beam 
and its absorption by decelerated atoms, the frequency 
detuning A*^j also essentially changes along the axis of 
the slower. 

The optimal offset of the equilibrium velocity v(z) from 
the resonant velocity is achieved at a point where the 
derivative of the force (I) reaches its maximum, because 
the damping of the relative motion of atoms around this 
point is maximal. It is easy to show that within our def- 
inition of the e coefficient (3) this optimal cooling condi- 
tion is achieved exactly at e = 0.75 at each local intensity 
of light inside the slower. The damping coefficient of the 
force, which is maximal at e = 0.75, reaches its maximum 
at so = 2 and slowly decreases with further increase of 
the intensity of the laser field. 

The calculation of the required spatial distribution of 
the magnetic field in a Zeeman slower proceeds as fol- 
lows. First, the actual velocity of the slowing atoms is 
calculated numerically according to the formula 
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With this approach it is easy to include the spatial 
variation of the saturation parameter sq(z) due to con- 
vergence of the cooling laser beam and even due to ab- 
sorption of its light by the decelerated atoms Q . Based 
on the calculated slowing velocity v z (z) , the correspond- 
ing distribution of the magnetic field can be derived from 
the equality A e //(z) = A*yj(z) as 



B(z) = h (S - kv(z) + A* eff (z)) /ii'. 



(0) 



This equation accepts multiple solutions, which corre- 
spond to different constant offcet of the magnetic field 
and corresponding frequency detuning Sq = A*^(zf) — 
kv(zf) — /j,'B(zf)/h of the slowing light, where Zf is the 
coordinate of the maximum of the magnetic field at the 
end of the slower. In this paper we are considering sym- 
metric spin-flip Zeeman slowers, in which the magnitudes 
of the magnetic field at both ends of the slowers are equal 
and of opposite sign. 
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III. LONGITUDINAL MD-ZEEMAN SLOWERS 
FOR SR ATOMS 

A. General approach to longitudinal MD-Zeeman 
slowers 

The field of a single magnetic dipole placed at the ori- 
gin of a Cartesian system of coordinates and oriented 
along the z-axis is described by the formulae 
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where r = \J x 2 + y 2 + z 2 and M is the magnetic moment 
of the dipole. 

In this paper the longitudinal Zeeman slowers based on 
arrays of magnetic dipoles, oriented parallel to the z-axis 
of the slower, are considered. 

The spatial distribution of the magnetic field gener- 
ated by an array of such dipoles is determined by sum- 
mation of the corresponding field amplitudes, produced 
by each of the individual dipoles. An array of mag- 
netic dipoles, oriented parallel to the z-axis and placed 
in the XOZ-plane with coordinates of individual dipoles 
of x = Xi, y — 0, z = Zi, are generating the magnetic field 
with distribution of the i? z -component at the z-axis of 
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FIG. 2: Spatial distribution of B z component of the magnetic 
field along the z-axis, generated by two continuous line-like 
arrays of magnetic dipoles, oriented parallel to that axis and 
located at the same distance R — 4 cm from it. The solid line 
corresponds to the short array, L s hort = 2 cm, and the dashed 
one to a long array of length Li ong = 40 cm. 
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FIG. 3: Simplest longitudinal spin-flip MD Zeeman slower for 
Sr atoms. 



For a symmetric distribution of the magnetic dipoles 
around the z-axis, all other components of the magnetic 
field at this axis are equal to zero. For positions away 
from the z-axis, the magnitudes of the transverse com- 
ponents of the field, B x and B y , can be computed in a 
similar way. 

The solid line in fig. 2 shows the spatial distribution 
of the B z (z) component of the magnetic field along the 
z-axis for a line-like continuous array of magnetic dipoles 
covering a length L s h or t = 2 cm with total magnetic mo- 
ment M = 3.1 A-m 2 , oriented parallel to this axis and 
placed at a distance of x = R = 4.0 cm from it. The 
dashed curve in fig. 2 shows the spatial distribution of 
the field at the z-axis generated by a Li ong= A0 cm long 
continuous line-like array of magnetic dipoles with total 
magnetic moment M — 62.0 A-m 2 , which is located at 
the same distance of x = R = 4.0 cm from the z-axis. 
It means, that the linear density of the magnetic mo- 
ments of these two arrays is the same. One can see that, 
in spite of the much larger total magnetic moment of 
the longer array of dipoles, its field amplitude is smaller 
than the field amplitude of the short array of magnetic 



dipoles. Moreover, for an infinitely long array of lon- 
gitudinal magnetic dipoles, the amplitude of its outside 
field becomes equal to zero. It is easy to show that the 
i? z -component of the magnetic field of such an array is 
maximal when its length is L ~ 0.4i?, while its spatial 
distribution along the z-axis is broadened but still has a 
single maximum. 

For laser cooling of Sr atoms, the 1 S'o — Pi transi- 
tion with wavelength A = 461 nm and natural lincwidth 
r = 2tt x 32 MHz is used, one of the simplest longitudinal 
spin-flip Zeeman slowers for Sr atoms, based on magnetic 
dipoles, is shown in fig. 3. It consists of four cylindrical 
neodymium magnets, placed symmetrically with respect 
to the z-axis of the slower in the X0Z plane. The mag- 
netic moments of the long and short magnets are directed 
parallel to the z-axis with opposite polarity. All the mag- 
nets have a diameter of 30 mm. The base of the two long 
magnets is located at z = and x = ±74 mm and their 
length is equal to 88 mm. The length of the short mag- 
nets is equal to 32 mm, while there distances from the 
axis of the slower are x = ±54 mm. The gap between 
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FIG. 4: The optimal magnetic field distribution (solid line) 
and the measured field of the transverse MD Zeeman slower. 



the magnets is equal to 65 mm. Therefore, the smaller 
end magnets are located closer to the axis of the slower, 
which provides a larger gradient of the magnetic field at 
the end of the slower. The total magnetic moment of 
each of the long magnets is Mi = B r V/[io = 54.45 A-m 2 , 
where V is the volume of the magnet, B r = 1.1 T is the 
residual induction of a standard neodymium magnet and 
fio = 4-7T • 10 _7 H-m _1 is the permeability of free space. 
The corresponding total magnetic moment of the short 
magnets is equal to = 19.8 A-m 2 . 

The central solid line in fig. 4 shows the target spatial 
distribution of the magnetic field, which was calculated 
for e = 0.9, capture velocity v c = 330 m/s, final velocity 
Vf = 25 m/s, slowing length of about 12 cm, a saturation 
parameter of the cooling laser beam at the end of the 
slower of sq — 1.0 and a convergence angle of the laser 
beam, 8 = 0.012. The design frequency detuning of the 
cooling laser beam is 5q/2tt = 434 MHz. The slower is 
designed for e = 0.9 to provide maximum deceleration 
of atoms, which minimises the length of the slower and 
decreases the transverse broadening of the atomic beam 
due to its heating by spontaneously scattered photons. 
It is worth mentioning that the operation of a Zeeman 
slower at almost critical deceleration (e = 0.9) is possible 
only when its spetial distribution of the magnetic field 
is smooth enough, as it is in the MD Zeeman slowers 
considered here. The thin solid line shows the magnetic 
field distribution along the axis of the slower, which was 
calculated with a MD model (8), where each magnet was 
substituted with a row of point-like magnetic dipoles, 
separated from each other by a distance of 4 mm, with 
the same total magnetic moment as the modeled mag- 
net. The lines are slightly shifted with respect to each 
other along the z-axis solely to see better the agreement 
between their slopes. The two dashed lines correspond 
to the spatial magnetic field distributions produced by 
each (left and right) half of the slower. One can see that 
the z-positions of maximum and minimum of the mag- 
netic field are located nearly exactly where the centres of 
the magnets are. Therefore, the slowing length of such a 
slower is shorter than the physical size of the slower and 
the magnets are not used efficiently. 
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FIG. 5: Design of the Sr longitudinal MD Zeeman slower. 
Each cylindrical magnet is 8 mm long and 20 mm in diameter. 



The design of a more efficient longitudinal Zeeman 
slower is shown in fig. 5. This slower consists of four 
identical sections, each of which includes two stacks of 
magnets ("entrance" and "output" magnets) of oppo- 
site polarity, placed symmetrically with respect to the 
z-axis and located in X0Z and Y0Z planes. Each "en- 
trance" magnet consists of 16 individual neodymium 
magnets with length of 8 mm, diameter of 20 mm and 
magnetic moment of M = 2.2 A-m 2 . Each "output" mag- 
net contains 6 individual magnets. The exact distances 
of the centres of these magnets from the z-axis, which 
are changing nearly linearly with z-coordinate, are given 
in Table 1. 

Unfortunately, the problem of finding of the positions 
of individual magnetic dipoles of the slower, which is pro- 
ducing the target spatial distribution of the field of the 
slower at its axis, has no single solution. The exact fitting 
of the target field distribution at the ends of the slower 
is also impossible. 

Our general procedure for finding the positions of the 
magnets consists of the following steps. First, a simple 
spin-flip slower with two long cylindrical magnets, placed 
parallel to the z-axis as in fig 3, is calculated to have its 
field distribution maximally close to the target field dis- 
tribution. As a next step, the magnets are moved closer 
to the axis of the slower, in such a way that the corre- 
sponding distances of the magnets Xi from the axis are 
linearly decreasing from the centre of the slower (where 
the magnetic field crosses zero) towards the outer ends 
of the slower. Finally, the residual kinks in the spatial 
distribution of the magnetic field are eliminated by fine 
tuning of the positions of the group of the magnets in the 
vicinity of that kink. 

This slower was designed for e = 0.92, capture veloc- 
ity v c = 410 m/s, final velocity vt — 25 m/s, satura- 
tion parameter at the end of the slower sq — 1.0 and 
diameter of the cooling laser beam at that location of 
1.0 cm with a convergence angle of 9 = 0.014 rad. The 
design frequency detuning of the cooling laser beam is 
6q/2tt — 476 MHz. The spatial distribution of the mag- 
netic field of the slower, calculated with the MD model, 
and the target field distribution are shown in fig 6. The 
slowing length is now 18 cm, which is close to the physical 
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FIG. 6: Spatial distribution of the magnetic field along the 
axis of the Sr longitudinal MD Zeeman slower. 
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FIG. 7: Spatial distribution of the magnetic field flux along 
the x-axis at the output of the Sr longitudinal MD Zeeman 
slower. 



length of the slower. 

To calculate the magnetic field distribution of the MD- 
Zeeman slower in the xy-plane one needs to take into ac- 
count all three components of the magnetic field B x , B y 
and B z , which define the total amplitude of the magnetic 

field \B(x,y,z)\ = 



/B% + B% + _B| responsible for the 

local Zeeman shift of atomic magnetic states. The trans- 
verse variation of the field is the strongest at the input 
and output edges of the Zeeman slower, where the dis- 
tances of the magnets from the axis of the slower are the 
smallest. The calculations show that within a region with 
diameter of 1 cm the field is changing by not more than 



TABLE I: Calculated positions of the magnets in the Sr 
slower. 
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1.5% (see fig. 7). The transverse component of the mag- 
netic field at the end of the slower and at x = ±0.5 cm 
and y = ±0.5 cm is only ±0.8 G. In the central region of 
the slower, where the distances of the magnets from the 
z-axes of the slower are larger, the transverse variation 
of the field is much smaller. 

In addition, a precise calculation of the MD-Zeeman 
slower based on finite-size neodymium magnets with di- 
ameter 20 mm and height 8 mm was performed with the 
finite element analysis software COMSOL. These calcu- 
lations show that the optimal distances of the centres of 
the magnets from the axis of the slower are exactly the 
same as in the MD model. 



IV. LONGITUDINAL MD-SLOWER FOR RB 
ATOMS 

Rubidium is a popular and convenient element for laser 
cooling, which is used in multiple atomic physics experi- 
ments, including producing BEC and laser-cooled atomic 
frequency standards. Therefore, there are many tradi- 
tional Zeeman slowers of different configurations that 
have been implemented in different labs. The Rb Zee- 
man slowers use the \5Si/ 2 , F — 2^ — > \bP 3 / 2 , F' = 3) 
optical transition, which has a wavelength of A = 780 nm 
and a linewidth, T/2tt = 6.07 MHz. The strength 
of this transition and the photon recoil momentum of 
the Rb atom are smaller than for Sr, which makes 
the typical length of a Rb Zeeman slower equal to 
1 m. Another problem of Zeeman slowing 87 Rb atoms 
is the crossing between the \5P a / 2 ,F' =3,m= —3) and 
\5Ps/2,F' = 2, m = — l) magnetic sublevels of the ex- 
cited 5P3/2 state, which takes place at B = 120 G. For the 
case of imperfect circular polarization of the cooling laser 



beam, the \5P 3/2 ,F' = 2 



T) state can be excited, 
which can decay into the "dark" ground hyperfine state 
\bS 1/2 ,F = l). The atoms in the \5S 1/2 ,F = l) state 
are out of resonance with the cooling light and therefore 
are lost from further deceleration. To preserve that from 
happening, Rb Zeeman slowers are usually is designed in 
such a way that they include an additional uniform offset 
field of about 200 G 

On the other hand, in the work [l2| the spin-flip 
longitudinal Zeeman slower has a field that changes 
from 200 G to — 100 G. This slower is using an addi- 
tional repumping light beam, which is tuned close to the 
|5S 1/2 ,.F = 1) \5P 1/2 ,F' = 2) transition on the D 1 
manifold and has the same frequency detuning as the 
main cooling laser beam. In this setup a record flux of 
slowed atoms of 3.2 x 10 12 at/s has been demonstrated. 

Our design of the spin-flip longitudinal Zeeman slower 
for 87 Rb atoms, based on permanent magnets, is shown 
in fig 8. The magnetic field of the slower is changing 
from -118 G to 116 G within its length of 0.56 m. Each of 
the magnets of the slower has a diameter 30 mm, height 
20 mm and magnetic moment M = 12.4 A-m 2 , while 
some of them are combined in longer cylindrical mag- 
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TABLE II: Calculated positions of the magnets in the Rb 
slower. 
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about 160 MHz. This makes it easy to generate the 
proper frequencies for the MOT and Zeeman slower cool- 
ing laser beams with standard single-pass acousto-optical 
modulators. 
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FIG. 9: Thick solid line - the target spatial distribution of the 
magnetic field of the Rb longitudinal Zeeman slower, calcu- 
lated for e = 0.8, capture velocity of v c — 275 m/s and final 
velocity Vf = 25 m/s. The dotted line shows the correspond- 
ing field distribution calculated with COMSOL software. The 
thin solid line is the field distribution calculated with MD 
model. 



V. CONCLUSION 



FIG. 8: Design of the Rb longitudinal MD Zeeman slower. 
Each cylindrical magnet is 20 mm long and 30 mm in diame- 
ter. 

nets. The exact positions of the centres of these magnets 
are given in Table 2. 

The slower is designed for e = 0.8, capture velocity 
v c = 275 m/s and final velocity Vf = 25 m/s. It is as- 
sumed that the cooling beam has a diameter at the end 
of the slower of 2.0 cm, where the saturation parameter 
is s = 3.0, and has a convergence angle of 9 = 0.02 rad. 
The design frequency detuning of the cooling laser beam 
is So/2tt — 196 MHz. The spatial distribution of the mag- 
netic field of the slower along its axis, calculated with the 
COMSOL software, is shown in fig 9. The sharp steps in 
the computed field are not real and are related to the er- 
rors due to the relatively large size of the mesh, which was 
used in the tree-dimensional finite-element computations. 
The last fact is confirmed by the MD-model simulations, 
which are in good agreement with the COMSOL results. 
The uniformity of the field in the transverse direction is 
very good. For example, at the output of the slower at 
z = 0.57m the increase in field at x = ±0.5 cm is about 
0.3 G only. 

One great advantage of the given spin-flip Zeeman 
slower for Rb is that its maximal Zeeman shift is only 



Examples of longitudinal spin-flip Zeeman slowers for 
Sr and Rb atoms, composed of arrays of compact perma- 
nent magnets, are presented. The advantage of these 
longitudinal slowers, compared to the transverse MD- 
Zeeman slowers, is that they demand only half as much 
power for the cooling laser light, which is important for 
some atoms, like Sr or Yb. In addition, the longitu- 
dinal spin-flip slowers work for all laser-cooled atoms, 
which is not the case for the transverse slowers. On 
the other hand, the longitudinal MD-Zceman slowers are 
more complicated, compared to the transverse ones, and 
demand larger magnets. The slowers do not consume 
electrical energy, and the spatial distribution of their 
magnetic field and the total length of the slowers are eas- 
ily adjustable. The calculated spatial distribution of the 
magnetic field of the slower is extremely smooth, which 
should make it possible to operate it at the optimal con- 
ditions, where the length of the slower is minimal, which 
reduces the transverse expansion of the slowed atoms. 



Acknowledgment 

Many thanks to Rachel Godun for the valuable com- 
ments. This work was funded by the UK National Mea- 
surement System Directorate of the Department of Trade 
and Industry. 



[1] W. Ketterle, Nature 435, (2005) 321. 



[2] E. W. Streed, A. P. Chikkatur, T. L. Gustavson, M. 



7 



Boyd, Y. Torii, D. Schneble, G. K. Campbell, D. E. 
Pritchard, and W. Ketterle, Rev. of Sci. Instr. 77, (2006) 
023106. 

[3] M. Takomoto, F. Hong, R. Higashi and H. Katori, Nature 

435, (2005) 321. 
[4] W. D. Phillips and H. Metcalf, Phys. Rev. Lett. 48, 

(1982) 596. 

[5] Y. B. Ovchinnikov, Opt. Comm. 276, (2007) 261. 
[6] Y. B. Ovchinnikov, Eur. Phys. J. Special Topics 163, 
(2008) 95. 



[7] Jabez McClelland, private communication. 
[8] P. Cheiney et al, ! arXiv:1101.3243l (20111. 
[9] I. Hill et al, to be published. 
[10] R. J. Napolitano, S. C. Zilio and V. S. Bagnato, Opt. 

Comm. 80, 110 (1990). 
[11] P. A. Molenaar, P. Van der Straten, H. G. M. Heideman 

and H. Metcalf, Phys. Rev. A 55, 605 (1997). 
[12] C. Slowe, L. Vernac and L. V. Hau, Rev. Sci. Instr. 76, 
(2005) 103101. 



